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ABSTRACT 


Interaction  of  a  small  system  S  with  a  large  reservoir  R  amounts  to  thermal 

relaxation  of  the  reduced  system  density  operator  P§(t).  The  presence  of  the 

reservoir  enters  the  equation  of  motion  for  Pg(t)  through  the  reservoir 

correlation  functions  f^-t)  (defined  in  the  text),  which  decay  to  zero  for  r  ■*  • 

on  a  time  scale  t^.  Commonly,  this  is  much  smaller  than  the  inverse 

relaxation  constants  for  the  time  evolution  of  Pg(t).  Then  a  series  of 

approximations  can  be  made,  which  lead  to  a  Markovian  equation  of  motion  for 

pg(t).  In  this  paper  the  assumption  of  a  small  reservoir  correlation  time  is 

removed.  The  equation  of  motion  for  pg(t)  i*  solved,  and  it  appears  that  the 

memory  effect,  due  to  £  0,  can  be  incorporated  in  a  frequency  dependence  of 

the  relaxation  operator  f(u).  Subsequently,  (unequal -time)  quantum  correlation 

functions  of  two  system  operators  are  considered,  where  explicit  expressions  for 

(the  Laplace  transform  of)  the  correlation  functions  are  obtained.  They  involve 

again  the  relaxation  operator  T(u),  which  accounts  for  the  time  regression. 

Additionally  it  is  found  that  an  initial-correlation  operator  f(w)  arises,  as  a 

consequence  of  the  fact  that  the  equal-time  correlation  functions  do  not 

factorize  as  Pg(t)  times  the  reservoir  density  operator.  It  is  pointed  out  that 

the  frequency  dependence  of  F(w)  and  the  occurrence  of  a  non-zero  f(w)  both  arise 

as  a  result  of  t  £  0,  and  should  therefore  be  treated  on  an  equal  footing. 

Explicit  evaluation  of  r(w)  and  f(w)  shows  that  their  matrix  eleswnts  can  be 

expressed  entirely  in  f^(w),  just  as  in  the  Markov  approximation.  Hence  no 

essential  complications  appear  if  one  should  go  beyond  the  limits  of  a  small 

reservoir  correlation  time  t  . 

c 

PACS:  02.90.4p,  05.40.4-J 
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I.  INTRODUCTION 

In  many  practical  casas  the  equation  for  the  evolution  of  the  density 
operator  p(t)  of  a  quantum  system  assumes  the  general  form 


i#dtp(t)  -  +  Hjj  +  Hj»p(t)]  . 


(1.1) 


where  H£  and  the  Hg  pertain  to  separated  components  S  (■  system)  and  R  (■ 
reservoir)  of  the  entire  configuration,  and  H£  denotes  an  interaction  between  S 
and  R.  Probably  the  most  familiar  example  is  spontaneous  decay  of  an  excited 
atom  in  empty  space.  Then,  H£  equals  the  atomic  Hamiltonian  (internal 
structure),  Hg  represent  the  electromagnetic  field,  and  HJ  is  the  dipole  coupling 
between  the  atom  and  the  electric  component  of  the  radiation  field,  which  causes 
the  spontaneous  transitions.  Since  Hg  has  a  large  (infinite)  number  of 
eigenstates,  an  exact  diagonalixation  of  the  complete  Hamiltonian  Hj  ♦  ♦  H|  ii 

intractable.  The  interest  is,  however,  in  the  behavior  of  the  atom,  as  it  is 
determined  by  its  interaction  with  the  radiation  field  (vacuum  or  black-body 
radiation).  Therefore,  one  introduces  the  reduced  atomic  (system)  density 
operator  by 


pg(t)  -  Trgp(t)  , 


(1.2) 


where  the  trace  runs  over  all  states  of  the  radiation  field  (the  reservoir).  The 
issue  of  reservoir,  relaxation  or  heat-bath  theory  is  then  to  derive  an  accurate 
equation  of  motion  for  Pg(t),  in  which  the  properties  of  R  only  enter  as  simple 
(and  explicit)  parameter  functions.  In  the  theory  of  spontaneous  decay  these  are 


the  Einstein  coefficients  and  the  Lamb  shifts. 
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Most  crucial  for  tha  development  of  a  relaxation  theory  is  the  concept  of  a 
large  reservoir.  If  the  system  S  were  not  present,  the  reservoir  would  be  in  a 
(thermal  equilibrium)  state  pR,  which  obeys 

[Hg.pg]  -  0  *  PR*  *  PR  .  TrRPR  *  1  »  (1.3) 

and  it  is  assumed  that  the  interaction  between  S  and  R  does  not  substantially 
affect  this  reservoir  state.  Or  more  precisely,  the  state  pR  changes  a  little 
due  to  the  interaction  with  S,  but  the  effect  on  the  time  evolution  of  the  system 
density  operator  Pg(t)  is  negligible.  In  the  quoted  example  this  implies  that  an 
atom  in  complete  vacuum  should  decay  in  the  same  fashion  as  an  atom  in  space  with 
a  single  photon  present.  As  a  consequence  of  this  large-reservoir  assumption,  we 
can  factorize  the  density  operator  as 

p(t)  z  PgtOPg  »  0-4) 

in  places  where  the  value  of  p(t)  determines  the  strength  of  the  interaction. 

In  order  to  derive  a  relatively  simple  equation  for  Pg(t),  a  sequence  of 

additional  approximations  is  usually  made,  which  rely  on  the  fact  that  the 

reservoir  correlation  time  t  is  short  in  comparison  with  the  inverse  relaxation 

c 

constants  1/r.  The  idea  is  as  follows.  One  derives  an  equation  for  P§(t)  which 
contains  a  quantity  of  the  form  <R(t)R(0)>,  with  R  a  typical  reservoir  operator 
(for  instance,  the  electric  field),  and  where  the  angle  brackets  indicate  an 
average  with  the  density  operator  pR,  e.g., 

<...>-  TrRpR(...)  .  (1.5) 
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Du*  to  th*  many  eigenvalues  Ku  of  and  the  large  cut-off  frequency  wc»  the 

function  <R(t)R(0)>  will  decay  to  zero  on  a  time  scale  of  the  order  of  t  -  1/w  . 

c  c 

On  the  other  hand,  as  a  result  of  the  interaction  between  S  and  R,  the  system 
density  operator  P§(t)  will  decay  on  a  time  scale  1/r  (with  r  an  Einstein 
coefficient,  for  instance),  and  in  many  cases  the  relation 


1*t  «  1  (1.6) 

c 

holds.  The  validity  of  (1.6)  allows  a  series  of  approx imat ions  (see  the 

Appendix),  commonly  referred  to  as  the  Markov  approximation. 

For  spontaneous  decay  the  restriction  (1.6)  is  rigorously  justified,  and  the 

equation  of  motion  for  Pg(t)  is  known  for  more  than  a  decade.1  ^  However,  not 

every  reservoir  has  a  short  correlation  time.  For  instance,  an  adsorbed  atom  or 

molecule  on  a  surface  interacts  with  the  substrate  through  phonon  coupling 

(crystals)  or  electron-hole  pair  creation  (metals).  In  the  case  of  physisorbed 

atoms  on  a  harmonic  crystal,  the  Hamiltonian  H£  accounts  for  the  kinetic  and 

potential  energy  of  the  atom.  The  potential  supports  bound  states,  separated  by 

>10  -  10  MHz  (infrared),  which  is  resonant  with  the  thermal  excitations  of  the 

crystal  (phonons).  Mechanical  coupling  (vibrations)  between  the  adsorbed  atom 

4-6 

and  the  lattice  atoms  gives  rise  to  thermal  relaxation  of  the  adbond  system. 
Typical  relaxation  constants  acquire  an  order  of  magnitude  of  10^  -  10^  MHz, 
whereas  the  cut-off  frequency  (Debye  frequency)  is  of  the  order  of  10*  MHz.  For 
electron-hole  pair  formation  the  situation  is  even  worse,  where  we  have  Ttc>1  so 
that  a  Markov  approximation  can  never  be  justified.^ 

There  exist  many  relaxation  theories.  Most  notable  are  the  projection 
techniques,^  ^  a  Langevin  formulation,^  and,  as  we  adopt  here,  a  reservoir 
approach.11  ^  A  feature  of  the  quoted  theories  is  that  they  all  lead  to  the 
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sane  result  as  soon  as  the  Markov  approximation  is  imposed.  Several  attempts 
have  been  undertaken  to  drop  this  Markov  assumption.*4  To  the  best  of  our 
knowledge,  however,  a  treatment  which  is  free  of  inconsistencies  and  leads  to 
explicit  expressions  (rather  than  formal  expressions  which  cannot  be  evaluated) 
was  never  formulated. 

II.  RESERVOIR  INTEGRAL 

In  this  section  we  set  up  the  notation  and  derive  an  integral  of  Eq.  (1.1), 
which  is  appropriate  for  imposing  the  reservoir  assumption.  The  first  step  is  a 
redefinition  of  the  system  Hamiltonian.  We  recall  that  the  interaction 
Hamiltonian  H^  is  an  operator  in  S+R-space,  and  therefore  its  reservoir  average 
<Hj>  will  be  an  operator  in  S-space.  In  order  to  eliminate  so-called  secular 
terms,  we  define  the  new  system  and  interaction  Hamiltonians  by 


Hs  -  HS  +  <Hi> 

• 

(2.1) 

Hj  •  Hj.  -  <Hj> 

• 

(2.2) 

and  the  advantage  of 

this 

rearrangement  comes 

from  the  fact  that  the  reservoir 

average  of  Hj  equals 

<ht>  »  0  . 

zero. 

Explicitly, 

(2.3) 

A  compact  and  transparent  representation  of  reservoir  theory  can  be  obtained 
with  a  Liouville-operator  formalism.  If  we  introduce  the  Liouvlllians  Lq  by 

L  o  -  if *(H  ,o]  ,  a  -  S,R, I  , 

a  a 


(2.4) 
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which  defines  their  action  on  an  arbitrary  operator  o  in  S+R-space,  then  the 
equation  of  motion  (1.1)  becomes 


i^p(t)  »  (Ls  +  1^  +  Lj)p(t) 


(2.5) 


Por  later  purposes  we  mention  a  few  properties  of  the  Liouvillians.  First,  Le 


and  commute,  since  they  act  on  a  different  part  of  Liouville  space.  Second, 
stands  for  a  commutator,  which  implies  the  relation 


Tr^o-d  , 


(2.6) 


for  any  o.  From  [H^Pp]  *  0,  Eq.  (1.3),  we  find 


e  PR  *  PR  * 


(2.7) 


and  due  to  the  shift  of  the  Interaction  over  its  average,  obeys 


Tr«Li(°sV  • 0  • 


(2.8) 


Bare  and  in  the  following,  Og  will  indicate  an  aribtrary  operator  in  S*space. 
An  integral  of  Eq.  (2.5)  reads 


-i(L_+LR)(t-t  ) 
p(t)  •  e  S  ^  °  p(tQ) 


•i: 


■i|  dt'  e 
t 

o 


i(Ls+LR)(t-f) 


LjP(t')  , 


(2.9) 


and  substitution  into  Eq.  (2.5)  then  yields 


.  -i(L_+LR)(t-t  ) 

i^p(t)  -  (Lg+L^)p(t)  +  Lre  S  K  °  P(tc) 

-ia.+LpKt-f) 

-iLjj  dt'  e  b  K  LlP(t')  ,  (2.10) 

*o 

which  is  an  exact  integral  of  the  equation  of  motion.  If  we  subsequently  take 
the  trace  over  the  reservoir  states,  the  left-hand  side  becomes  idpg(t)/dt,  which 
equals  the  rate  of  change  of  the  system  density  operator  due  to  the  free 
evolution  (the  term  LgPg(t)  on  the  right-hand  side)  and  the  coupling  to  the 
reservoir  (terms  proportional  to  L^).  Hence  the  integral  in  Eq.  (2.10)  accounts 
for  the  relaxation  of  Pg(t),  and  its  value  is  proportional  to  the  coupling 
strength.  Therefore  we  can  adopt  the  reservoir  assumption,  Eq.  (1.4),  on  p(t') 
in  the  integrand.  We  then  find  the  equation  of  motion  for  Pg(t)  to  be 

.  -i(L-+LB)(t-t  ) 

idtps(t)  “  LSpS(t)  +  TrRLI*  p(to} 

ft  -i(Lc+Lp)(t-t' ) 

-iTr^Ljl  dt'  e  L^PgU'^)  ,  (2.11) 

*o 

for  t  >  tQ.  It  is  important  to  note  that  the  initial  value  p(tQ)  of  the  density 
operator  (not  the  system  part)  remains  present  in  the  equation  of  motion  for 
pg(t),  in  general.  Equation  (2.5)  determines  the  time  evolution  of  p(t)  for 
t  >  tQ,  and  the  solution  of  Eq.  (2.5)  is  fixed  as  soon  as  an  initial  value  p(tQ) 


I 

\ 


is  prescribed.  Since  p(t  )  is  not  determined  by  the  equation  of  motion,  a 

o 

further  specification  of  the  initial  state  p(tQ)  is  necessary. 

III.  DENSITY  OPERATOR 

For  finite  memory-time  reservoirs  the  choice  of  p(tQ)  is  more  than  a  matter 

of  convenience.  If  the  system  has  been  in  contact  with  the  reservoir  prior  to 

t  ,  then  p(t  )  is  determined  by  its  time  evolution  in  the  recent  past  t  <  t  ,  and 
o  o  o 

consequently  the  value  of  p(tQ)  is  no  longer  arbitrary.  As  a  solution,  we  simply 

define  the  instant  of  time  t  as  the  time  point  at  which  the  interaction  Lr  is 

O  X 

switched  on.  We  can  then  always  take  t  to  be  arbitrarily  far  into  the  past. 

For  t  <  t  the  reservoir  is  in  its  thermal-equilibrium  state  pD,  and  the  system 

density  operator  Pg(t)  evolves  independently  of  the  reservoir.  Therefore,  we 

have  for  t  <  t 
-  o 

p(t)  *  pg(t)pR  .  (3.1) 

Substitution  into  Eq.  (2.11)  and  applying  Eqs.  (2.7)  and  (2.8)  then  shows  that 
the  term  with  p(tQ)  vanishes  identically,  due  to  the  shift  of  the  interaction 
Hamiltonian  over  its  average.  Then  the  equation  of  motion  for  Pg(t)  becomes 

^dt^S^^  "  LSPS^ 

ft  -KL.+LaKt-t1) 

-iTr^jj  dt'  e  K  L^p^f)^)  ,  (3.2) 

*o 

for  t  >  t  . 

-  o 

Solving  Eq.  (3.2)  is  most  easily  done  in  the  Laplace  domain.  If  we  define 
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(w)  ■  f 

Jt 


iw(t-t  ) 


ps(t)  . 


(3.3) 


then  the  transformed  equation  of  motion  reads 


(»-Ls)5s(»)  -  ips(t0)  -  ITrgLj  . 


(3.4) 


with  solution 


PS(*>  "  w-Ls+ir(u))  pS^o^ 


(3.5) 


Here  we  introduced  the  relaxation  operator  f(w)  as 


f(“)0S  -  TtRLI  V°sV  • 


(3.6) 


which  can  equivalently  be  written  as 


i(w-L--L_h 

f(w)°S  "  TrRLlJ  dT  *  LI(osV 

*  o 


(3.7) 


From  Eq.  (3.5)  we  see  that  p c(w),  and  thereby  p_(t)  for  t  >  t  ,  is  determined  by 

a  a  —  o 

Pc(t  )  only,  and  not  by  p_(t)  for  t  <  t  .  This  is  of  course  a  result  of 

O  O  j  O 

assumption  (3.1).  The  memory  in  the  time  evolution  of  Pg(t)  is  displayed  in  the 
frequency  dependence  of  f(w).  In  the  Appendix  we  show  that  F(io)  acquires  a 
constant  value  (w- independent)  in  the  Markov  approximation. 

From  Eq.  (3.7)  we  notice  that  f(cj)  has  the  form  of  a  Laplace  transform 
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F(u>)  -  dT  •1WT  I*(t)  , 


(3.8) 


where  T(r)  Is  given  by 


r(t)os  -  TrRLj  e'i(LS+LR)x  Lj(<>gPR)  , 


(3.9) 


for  t  >  0.  Rewriting  the  equation  of  motion  (3.2)  in  terms  of  T(t)  gives 


idtpS^t^  *  LSPS^ 


-if  dt' 

Jt 


r(t-t')ps(t')  , 


(3.10) 


which  reveals  that  the  time  width  of  T(t)  (its  decay  time  for  t  >  0)  equals  the 
Demory  time  of  the  reservoir- interaction  term.  It  is  the  width  of  r(t)  which  is 
usually  termed  the  reservoir  correlation  time  t  .  Then  it  follows  from  Eq.  (3.8) 
that  the  frequency  width  of  F(w)  is  of  the  order  of  1/tc,  and  for  -*■  0,  F(u) 
becomes  independent  of  w. 

The  time  evolution  of  p_(t)  for  t  >  t  will  have  little  significance  in 

w  ""  O 

general,  which  is  partially  due  to  the  factorization  at  t  *  tQ.  Due  to  the 
coupling  to  the  reservoir,  the  density  operator  Pg(t)  will  relax  to  a  steady 
state  (thermal  equilibrium) 


Ps  ■  lim  Pg ( t ) 
t-*» 


(3.11) 


on  a  time  scale  1/r,  as  mentioned  in  the  Introduction.  Here,  T  denotes  a  typical 
matrix  element  of  f(w)  (not  of  T(t)).  From  the  identity 
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p  -  lim  -lu)ps(u.)  (3.12) 

u r+0 

and  Eq.  (3.5),  we  find  the  equation  for  pg  to  be 

(Ls-iP(0))ps  «  0  .  (3.13) 

This  shows  that  the  long-time  solution  of  pg(t)  is  determined  by  f(u>)  at  u  *  0. 
Furthermore,  we  notice  that  the  dependence  on  the  initial  value  p„(t  )  has 

w  O 

disappeared  in  Eq.  (3.13),  which  reflects  that  the  memory  of  the  preparation  of 
the  system  at  t  is  erased. 


IV.  CORRELATION  FUNCTION 

Measurement  of  the  steady-state  density  operator  pg  of  a  physical  system  is 
tantamount  to  the  determination  of  its  relaxation  constants,  which  are  the  matrix 
elements  of  f(o),  as  displayed  in  Eq.  (3,13).  Dynamical  properties  of  the  system 
in  contact  with  the  reservoir,  however,  are  reflected  in  the  time  evolution  of 


Pg(t)  before  it  reaches  its  steady  state  Pg.  In  view  of  Eq.  (3.5),  this 


transient  behavior  of  Pg(t)  is  governed  by  the  frequency  dependence  of  the 
relaxation  operator  F(w).  Besides  the  fact  that  a  density  operator  is  not 
amenable  to  direct  observation,  we  also  see  from  Eq.  (3.5)  that  the  details  of 
pc(t)  depend  on  the  preparation  of  the  system  at  t  ■  t  .  Obviously,  it  is 

2)  O 


impossible  to  fix  pc(t  )  (say,  the  wave  function  of  an  atom)  at  a  single  Instant 

w  O 


of  time,  and  subsequently  measure  its  evolution  for  t  >  t^. 


A  standard  method  of  obtaining  dynamical  information  about  a  system  is  by 
observation  of  steady-state  correlation  functions  of  system  operators,  say  X  and 


Y.  If  we  take  arbitrarily  t  as  the  instant  of  time  at  which  the  Schrijdinger  and 


Heisenberg  pictures  coincide,  then  the  time  dependence  of  X  is  given  by 


. .  ~  y -  /■/-  l  -  •  -  ’  > .  -  --  -  •  ^  i,  *•' ■  ~  - 


•.<  .1  ■<»  .1  ’.I  *.»  J. 


I.  ■  J  .!»  J-.J-lJI  . 


x(t)  -  e 


iL(t-t  ) 


(4.1) 


where  L  Indicates  the  Liouvillian  Lg+L^+Lj  of  the  entire  system.  Then  X(tQ)  *  X 
is  an  operator  in  S*space  only,  but  for  t  >  tQ,  X(t)  is  an  operator  in  S+R-space, 
due  to  Lj.  Hence  the  time  evolution  of  X(t)  carries  information  on  the 
interaction  with  the  reservoir.  The  correlation  function  of  two  operators  X  and 
Y  is  defined  as  the  expectation  value 


«X(t')Y(t)»  -  Trp(t  )X(f)Y(t) 

o 


(4.2) 


The  double-bracket  notation  indicates  an  average  with  the  full  density  operator 
of  S+R,  rather  than  with  only  (Eq.  (1.5)).  Transformation  of  Eq.  (4.2)  to  the 
Schrodinger  picture  gives 


«X(t')Y(t)»  -  TrY  e‘iL(t"t')(p(t,)X)  , 


(4.3) 


or  equivalently 


«X(t')Y(t)»  -  TrX  e-1L<t,-t)(Yp(t)) 


(4.4) 


We  notice  that  the  initial  time  tQ  has  disappeared  in  Eqs.  (4.3)  and  (4.4),  which 
already  removes  the  ambiguities  associated  with  the  preparation  of  p_(t  ).  A 
steady-state  correlation  function  is  now  defined  as  <<X(t')Y(t)»  with  t  >>  tQ, 
t*  >>  t  and  t-t'  fixed.  Then  the  system  is  in  state  pc,  which  is  time 
independent  and  a  solution  of  Eq.  (3.13).  The  time  regression  of  the  correlation 
functions  (their  t-t' -dependence)  is  governed  by  the  same  exponential  which 
determines  the  time  evolution  of  p(t),  and  therefore  we  can  extract  dynamical 


r 
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properties  of  the  system  by  an  observation  of  the  steady-state  correlation 
functions. 

Commonly,  time  regressions  are  not  measured  directly.  For  atoms  or 
molecules  on  a  solid  substrate,  for  instance,  one  determines  the  spectral  profile 
for  the  absorption  of  low-intensity  monochromatic  laser  radiation  with  frequency 
u>.  The  spectral  distribution  as  a  function  of  u  and  in  the  steady  state  is  then 
given  by  expressions  of  the  form 

I(w)  -  lim  f"  dt'  eiw(t''t)«X(t,)Y(t)»  ,  (4.5) 

t-H»  Jt 

which  will  further  be  referred  to  as  the  spectrum.  It  is  the  goal  of  this  paper 
to  evaluate  I(w)  for  a  system  in  interaction  with  a  finite  memory-time  reservoir. 

V.  SPECTRUM 

From  Eq.  (4.5)  we  observe  that  we  need  the  correlation  function  for  t*  >  t, 
and  therefore  the  representation  (4.4)  is  most  suitable.  Then  the  occurring 
exponential  is  the  same  as  for  the  time  evolution  of  p(t).  If  we  introduce  the 
Hilbert-space  operator  (Liouville-space  vector) 

A(t'.t)  -  e'iL(t''t)(Yp(t))  ,  (5.1) 

then  the  correlation  function  can  be  represented  by 

«X(t '  )Y(t)»  -  TrsXAs(t\t)  , 


(5.2) 


u 


which  only  involves  the  systea  pert  A^t'.t)  •  TrRA(t',t).  In  tens  of  the 
Laplace  transfon  with  respect  to  t* 

As(w,t)  -  £df  elta,(t'’t)As(t,,t)  ,  (5.3) 

the  spectrum  attains  the  fon 

I(w)  -  lim  Tr_XAc(«n,t)  .  (5.4) 

t—  5  5 

Differentiating  Eq.  (5.1)  with  respect  to  t*  yields  the  equation  of  motion 
for  A(t' ,t) 

ift7  A(t'.t)  -  (Ls  +  Lj  ♦  LJ)A(t',t)  ,  (5.5) 

which  has  to  be  solved  for  t'  >  t,  with  initial  value 

A(t.t)  •  Yp(t)  .  (5.6) 

Equation  (5.5)  is  identical  to  Eq.  (2.5)  for  p(t),  and  integrals  can  be  found  in 
the  same  way.  The  difference  between  a  density  operator  and  a  correlation 
function  ls  that  for  p(t)  we  can  choose  the  Initial  value  p(tQ)  arbitrarily, 
whereas  for  A(t',t)  the  initial  value  ls  unambiguously  given  by  Eq.  (5.6).  This 
reflects  the  fact  that  A(t',t)  is  essentially  a  two-time  quantity.  Its 
regression  from  t  to  t'  is  governed  by  Eq.  (5.5)  and  its  dependence  on  t  enters 
through  the  initial  condition,  Eq.  (5.6).  The  memory  in  the  time  regression,  due 
to  the  finite  reservoir  correlation  time,  is  of  course  the  same  as  for  the 


IS 


density  operator  and  can  be  accounted  for  by  the  frequency*dependent  relaxation 
operator  f(u),  as  we  shall  show  below.  As  a  second  effect  of  a  finite  tc  the 
density  operator  p(t)  In  the  Initial  value  will  carry  a  memory  of  Its  time 
evolution  in  the  recent  past.  It  is  tempting  to  argue  that  we  consider  the 
steady  state  t  -*  •,  so  that  the  density  operator  p(t)  is  constant  in  tine.  By 
the  large-reservoir  assumption  we  know  that  the  reservoir  remains  in  the  state 
pg,  whereas  the  system  is  in  state  pg  for  t  ■»  •.  This  would  imply  the 
replacement  p(t)  -*  p  pg  in  Eq.  (5.6),  which  in  turn  would  eliminate  the  explicit 
t-dependence  of  A(t',t),  making  the  limit  t  -*  •  in  Eq.  (5.4)  trivial.  We  shall 
show  that  this  procedure  cannot  be  justified  if  is  finite. 

Since  Eq.  (5.5)  is  identical  to  Eq.  (2.5)  for  p«.(t),  we  can  derive  the 
appropriate  integral  along  the  same  lines.  The  analogue  of  Eq.  (2.11)  is 

.  -i(l  +L_)(t'-t) 

i^T  As(t\t)  -  LsAs(f,t)  ♦  TrRLj  e  *  (Yp(t)) 

,t'  -i(L_+Lp)(t' -t") 

-iTrRLJ  dt"  e  5  K  LI(As(t",t)pR)  .  (5.7) 


which  contains  p(t)  explicitly.  Now  we  can  substitute  the  right-hand  side  of  Eq. 

(2.9)  for  p(t)  and  take  for  p(t  )  the  value  p_(t  )p_.  Then  Eq.  (5.7)  becomes 

o  a  o  K 


ldt'  "  LsAs^t' 


+  TrRLj  e 


-i(Lg+LR)(t' -t)  -i(Ls+LR)(t-tQ) 

J  h  • 


(pS(to,5*) 


■iTr, 


-i(L_+LB)(t'-t)  ft  -i(L  +Lj(t-t") 
e  b  K  Ly  j  dt"  e  5  K  LlP(t") 


16 


ft'  -i(L-«V<f-t") 

•ITrjjLj  j  df  •  5  *  LI(As(tM,t)pR)  .  (5.6) 


where  we  introduced  tha  Liouvilllan  Ly  by 


Vs  "  YoS  • 


(5.9) 


in  ordar  to  avoid  notations  with  too  many  brackats.  In  tha  sacond  tars  on  tha 
right-hand  sida  of  Eq.  (5.8),  the  exponentials  with  Lg  act  only  on  pR,  because  Lg 
coamutes  with  and  Ly.  Therefore  they  cancel,  according  to  Eq.  (2.7).  The 
renaining  two  exponentials  and  Ly  affect  only  ps(tQ),  and  the  result  is  soaa 
operator  og  in  S-space.  With  Eq.  (2.8)  we  then  find  that  the  whole  term  is 
identically  zero.  Considering  the  third  tarn  on  the  right-hand  side,  we  notice 
that  it  has  the  form  of  a  reservoir  integral,  as  in  Eq.  (2.10),  which  implies 
that  we  can  factorize  p(t")  here.  Then  we  define  a  'density  operator'  pSR(t)  of 
S+R- space  by 


-i(Lc+LB)(t-t') 

{(t)  »  J  dt'  c  Lj(ps(t')pR)  , 

t.\ 


(5.10) 


which  allows  us  to  write  Eq.  (5.8)  as 


idt'  "  LSAS^t'*t^ 


"iTrpLj  • 


-i(Ls+LR)(t'-t) 


VsR(t) 


,t'  -  i(L_+L,j)(t ' -t") 

jLjJ  dt"  e  b  *  LI(As(t",t')pR)  .  (5.11) 


rvW 


A  *-  %  ,N  *.  . 
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Kext  v*  take  the  Laplace  transfora  of  Iq.  (5.11),  recalling  that 


Ag(t.t)  •  LyPs(t) 


(5.12) 


as  follows  froa  Eq.  (5.6),  and  rearrange  the  teras.  We  then  obtain 


Vw,t)  "  w-Ls+ir(w){LYpS(t)  '  Tr»LI  w-Lj-Lj  *  (5-13) 


The  factor  in  front  of  the  curly  brackets  is  the  saae  as  in  Eq.  (3.5),  and  it 


represent  the  tiae  regression  from  t  to  t*  of  As(t',t).  The  first  tera  inside 


the  brackets,  L^PgCt),  corresponds  to  a  factorised  initial  state.  If  we  would 
have  replaced  A(t,t)  ■  Yp(t)  by  Y(ps(t)pR),  then  it  is  easy  to  see  that  the 
second  tera  on  the  right-hand  side  of  Eq.  (5.7)  would  have  disappeared,  and 
thereby  the  second  tera  in  curly  brackets  in  Eq.  (5.13).  Conversely,  the  tera 
with  PSR(t)  in  Eq.  (5.13)  accounts  for  the  correlations  between  S  and  R  in  p(t), 
which  are  present  at  the  initial  tiae  for  the  tiae  evolution  of  A(t',t)  froa  t  to 
t'  . 

The  explicit  tiae  dependence  of  l.(w,t)  enters  through  p  (t)  and  PCD(t).  If 

^  5  SR 

we  denote  their  steady-state  values  by  an  overbar,  then  the  spectrua,  Eq.  (5.4), 

becoaes 


K“)  "  ^sbe  w-Lj'f il*(w)  {LYPS  ’  TrRLI  w-Lg-LRLYpSl) 


(5.14) 


with 


Xo, 


(5.15) 


A 


A 


A 
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Expression  (5.14)  Involves  the  systea- reservoir  stste  psr,  which  night  seen 

cuabersosM.  Froa  Eq.  (5.10)  we  find  the  Leplece  trsnsforn  of  pgg(t)  to  be  j 

( 
< 

(5.i6)  : 


in  terns  of  Ps(u)  froa  Eq.  (3.5).  Then  the  steady- state  psr  follows  froa  the 
identity  (3.12),  which  gives 


iSR<“)  ‘ 


s*  «S-hi 


VW 


(5.17) 


Here,  the  notation  10+  Indicates  a  small  positive  imaginary  part,  which  is 
necessary  to  assure  the  convergence  of  Laplace-transfora  integrals,  or 
equivalently,  the  existence  of  the  Inverse  of  10+-LS 
show  how  to  evaluate  the  right-hand  side  of  Eq.  (5.17).  If  we  define  an  operator 
T(w)  by 


Lg.  In  the  next  section  we 


tU)os  •  Tr„L 


*“1  -L_  Sf 


10  -1.3-4 


Li  (Vf) 


(5.18) 


then  the  spectrua  attains  the  form 

I<“>  -  T'sS  <V‘f<"))5s  •  <5  t,) 

Equation  (5.19)  is  the  most  condensed  and  general  representation  of  the  result  of 
this  paper.  The  finite  memory-time  of  the  reservoir  appears  as  a  frequency 
dependence  of  the  relaxation  operator  T(u>),  and  as  a  non-vanishing  initial- 
correlation  operator  T(u>). 
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VI.  INTERACTION 

Although  the  result  (5.19)  is  appealing  and  explicit,  the  occuring  operators 
T(u)  and  f(u>),  which  represent  the  interaction  of  the  systea  with  the  reservoir, 
night  look  awkward  in  their  definitions,  Eqs.  (3.6)  and  (5.18).  Especially  the 
reservoir  Liouvi Ilian  in  denominators  and  the  appearance  of  i0+  in  £q.  (5.18) 
night  seen  to  make  an  explicit  evaluation  of  f(w)  and  ?(u)  intractable.  Such  is 
however  not  the  case,  as  we  shall  show  in  this  section. 

Obviously,  an  elaboration  of  f(u>)  and  T(u>)  requires  additional 
specifications  of  the  interaction  Hamiltonian  Hj.  It  will  turn  out  to  be 
sufficient  to  assume  the  form 

Hj  -  A  £  £k#k  .  (6.1) 

k 

with  Sk  <<^>  a  pure  S-  (R-)operator.  The  form  (6.1)  pertains  to  most  practical 
situations  we  have  encountered.  In  the  case  of  fluorescence,  signifies  the 
k~th  Cartesian  component  of  the  atomic  dipole  moment,  and  for  adsorbates  on  a 
substrate  the  subscript  k  takes  on  two  values,  corresponding  to  the  two  terms  in 
the  binding  (Horse)  potential.  In  fact,  the  form  (6.1)  for  Hj  can  always  be 
enforced  by  an  expansion  in  matrix  elements. 

Evaluation  of  the  relaxation  operator  f(u)  starts  from  its  representation 
(3.9)  in  the  time  domain.  We  expand  the  two  Lj's  **  commutators,  which  gives 
rise  to  four  terms.  Then  we  insert  Hj  from  Eq.  (6.1),  and  we  notice  that  every 
factor  is  an  operator  in  S-  or  R-space  only.  Combining  the  R-operators  and 
taking  the  trace  over  the  reservoir  states  then  shows  that  the  R  contribution  can 
be  accounted  for  by  a  single  complex-valued  function 
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fkl(o  •  <$k«  R  tHt>  • 


(6.2) 


which  will  be  celled  the  reservoir  correlation  function.  We  find 


-1LsT 


r(,)oS*ILl*  S  •  f*k(T)°s4>  - 


(6.3) 


kt 


with 


Lt°s  "  * 


(6.4) 


Expression  (6.3)  only  involves  the  system  operators  and  the  Liouvillian  Lg 

% 

for  the  free  evolution  of  the  system.  The  reservoir  enters  via  the  parameter 


functions  fRt(x),  which  can  be  found  as  soon  as  a  particular  reservoir  is 


prescribed.  For  a  harmonic  crystal,  for  instance,  the  reservoir  correlation 
functions  are  given  analytically  in  Ref.  17. 

The  initial  correlation  operator  ?(u)  from  Eq.  (5.18)  is  the  Laplace 
transform  of 


T(t)os  -  TrRLj  • 


'1(LS+LR)t.  1  .  ,  -  . 

hr  _+  .  .  li(ospr)  • 


(6.5) 


10  -S'hi 


First  we  recall  that  the  notation  iO  should  be  read  as 


iO 


1  -  C  i(w“Lg-L_)x ' 

7— —  li(V«)  '  **  u*  J  d'’  •  li(Vr> 


(6.6) 


ur»i0 


2 


Then  we  insert  the  form  (6.1)  for  the  interaction  and  rearrange  the  S-  and  li¬ 
terals.  We  then  obtain  for  T(t) 


t(t)os 


dt'  e 


■V 


« (v,+,,,^kos  •  • 


(6.7) 


where  the  reservoir  is  again  entirely  incorporated  in  the  functions 

Since  the  f^tO's  decay  to  zero  sufficiently  fast  for  t  ■*  »,  we  omitted  at  this 

stage  the  i0+  in  the  right-most  exponential. 


VII.  LAPLACE  TRANSFORM 

Before  we  can  take  the  Laplce  transform  of  Eqs.  (6.3)  and  (6.7),  we  must 
work  out  the  exponentials.  Eigenstates  of  the  system  Hamiltonian  will  be  denoted 
by  |a,a>,  where  a  indicates  the  energy  and  a  any  degeneracy.  By  definition  they 
obey 

Hg |a,a>  «  hw^| a,a>  (7.1) 

With  respect  to  its  own  eigenstates  we  can  write  Hg  as 

Hg  -  £  h(»a|a,a><a,a|  -  £  »  (7.2) 

a.a  a 


with 
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P#  *  £  |a,a><a,aj  , 


(7.3) 


the  projector  on  the  subspace  with  energy  #u>a.  From  the  orthonoraality  of  the 
states  |a,o>  we  have 


P  P.  -  6  .  P  , 
a  b  ab  a 


(7.4) 


and  from  the  completeness  of  the  set  ja,a>  we  find  the  closure  relation 


(7.5) 


Then  it  is  an  easy  matter  to  expand  the  exponential  exp(-iLgT)  in  projectors, 
which  gives 


-1LsT 

e  o. 


r  . t 

b 


°spb  • 


(7.6) 


in  terms  of  the  level  separations 


4  .  *  w  -  <u. 

ab  a  b 


(7.7) 


Next  we  subsitute  Eq.  (7.6)  into  Eq.  (6.3)  and  evaluate  the  Laplace 
transform.  We  obtain 


r<“)os  •  l  Lt  l  p.(ftk(4b.+“,*V’s  -  F*k(-4b.-"),’s*Vpb  • 


kt  ab 


(7.8) 
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in  terns  of  the  Laplace  transform  of  f ^ t ) .  Because  equals  the 

coamutator  with  eP^,  the  right-hand  side  only  involves  operators  and 
projectors.  If  we  insert  the  closure  relation  (7.5)  in  various  places  in  Eq. 
(7.8),  we  immediately  find  the  matrix  representation  of  f(w)  in  terns  of  matrix 
elements  of  The  result  (7.8)  is  the  most  compact  representation  of  the 

explicit  form  of  f(ui). 

In  the  very  same  way  we  find  the  Laplace  transform  of  T(t)  from  Eq.  (6.7), 
although  with  considerably  more  effort,  which  is  due  to  the  double  integral  (over 
t  and  t').  The  result  is 


T  (u))os 


kt  abc 


1 


A  +uj 
ac 


P  YP 
c  a 


‘  (fJk('Ab,)  '  !tk<'4bc'“))oS^k,Pb 


(7.9) 


which  has  a  striking  resemblance  with  Eq.  (7.8).  Most  remarkable  is  that  T(u>) 
can  again  be  expressed  in  the  reservoir  correlation  function  fjct(UJ)  which  also 
determines  F(w),  and,  as  shown  in  the  Appendix,  the  relaxation  operator  in  the 
Markov  approximation.  The  distinction  is  that  the  functions  £^(01)  occur  with 
different  arguments. 

If  w  equals  a  level  separation  A£a  •  -A  ,  then  the  denominator  of  the  first 
factor  under  the  triple  summation  in  Eq.  (7.9)  becomes  zero.  For  u>  *  Aca  we  have 
Abc  +  w  *  A^,  an<*  hence  the  difference  of  the  two  functions  J  ^  in  curly 
brackets  also  approaches  zero.  In  the  process  of  deriving  Eq.  (7.9)  we  found 
that  this  feature  does  not  constitute  a  problem.  The  limit  is  simply 
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1 


limA  4  +w'*tk'“ba'  *tk'“bc' 

w*-A  ac 
ac 


SI  ftk<">  •  ln  “  ■  V  • 


(7.10) 


and  there  is  no  singularity  or  discontinuity  if  w  passes  across  a  resonance. 


VIII.  CONCLUSIONS 

If  the  decay  time  x  of  the  reservoir  correlation  function  f..(x)  for  x  ■*  • 

C  Hi 

is  not  small  in  comparison  with  the  relaxation  times  1/T,  which  are  determined  by 

the  same  function  (see  Eq.  (7.8)),  then  a  Markov  approximation  cannot  be  correct. 

In  this  paper  we  imposed  no  limits  on  X£.  We  only  assumed  that  the  system  S  is 

snail  in  comparison  with  the  reservoir  R.  The  finite  value  of  X£  amounts  to  a 

memory  in  the  time  evolution  of  the  density  operator,  which  is  reflected  in  a 

frequency  dependence  of  the  relaxation  operator  f(w).  Correlation  functions  of 

system  operators  depend  on  two  times,  t’  and  t.  The  regression  from  t  to  t' 

exhibits  the  same  memory  effect  as  the  time  evolution  of  the  density  operator. 

Additionally,  the  equal-time  correlation  function,  which  is  the  initial  value  for 

the  time  regression,  carries  a  memory  to  the  recent  past.  It  appears  that  this 

second  phenomenon  could  be  accounted  for  by  an  initial  correlation  operator  T(u>) 

in  the  expression  for  the  spectrum  I(w). 

Frequency-dependent  relaxation  operators  are  widely  applied  in  the 

literature.  Their  Laplace  inverse  r(x)  is  sometimes  called  a  memory  kernel, 

because  it  is  the  finite  time-width  of  T(x)  which  brings  about  the  memory  in  the 

time  evolution,  as  is  most  obvious  from  Eq.  (3.10).  Initial-correlation 

18 

operators,  however,  are  rare.  Despite  the  fact  that  the  frequency  dependence 
of  F(uj)  originates  from  the  same  memory  mechanism  which  amounts  to  a  non¬ 
vanishing  T(u>),  the  latter  is  usually  not  found.  As  pointed  out  in  the 


derivation  of  f(w),  the  disappearance  of  T(u)  is  a  consequence  of  a  factorization 
of  the  initial  value  or  state,  which  cannot  be  justified  in  general. 
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correlation  time  of  a  stochastic  process,  which  drives  a  multiplicative 
stochastic  differential  equation. 
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APPENDIX.  MARKOV  APPROXIMATION 

It  is  illuminating  to  compare  the  results  of  the  present  paper  with  their 

equivalents  under  the  Markov  approximation.  Since  we  already  have  Eq.  (2.11)  we 

can  start  there.  First,  we  state  that  we  are  not  interested  in  a  time  evolution 

of  pc(t)  on  a  time  scale  t  ,  which  implies  that  we  can  factorize  p(t  )  as 
3  c  o 

P«(t0)pjj.  In  this  fashion  we  discard  the  memory  of  the  initial  state  to  its 

past,  which  immediately  gives  f(w)  *  0,  or  equivalently,  the  second  term  on  the 

right-hand  side  of  Eq.  (2.11)  is  zero.  Second,  we  know  that  if  we  work  out  the 

integral  in  Eq.  (2.11),  we  find  reservoir  correlation  functions  f^(t-t'),  which 

decay  to  zero  on  a  time  scale  t  .  Therefore,  the  major  contribution  to  the 

c 

integral  comes  from  t  -t  <  t'  <  t.  Because  we  impose  the  condition  T t  <<  1, 

c  -  —  c 

the  density  operator  Pg(t')  in  the  integrand  is  not  affected  significantly  by  the 
relaxation  process  on  this  small  time  interval.  Then  we  can  replace  p^t')  by 
its  free  evolution 

iL  ( t- 1 '  ) 

Ps(t' )  ■  e  Pg(t)  ,  (A1 ) 

and  subsequently  take  p,.(t)  outside  the  integral.  Third,  according  to  the  first 
assumption  we  can  take  t-t  >>  which  gives  in  combination  with  the  fact  that 
the  integrand  is  only  non-zero  on  a  time  interval  tc  that  we  can  replace  tQ  by 
minus  infinity.  Combining  everything  then  yields 

idtPS^  "  ^LS*irM^PS^  ’ 


with 
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-‘W' ,  ‘V,  - . 

*  LI  *  °S°R^ 


The  Laplace  transform  of  Eq.  (A3)  reads 


rM°S 


TrRL 


I  f dT 

*  A 


ps(ta,) 


w-LS+irM 


ps(V 


» 


(A3) 


(A4) 


and  comparison  with  Eq.  (3.S)  then  shows  that  is  the  Markovian  equivalent  of 

T(ui),  and  indeed,  the  frequency  dependence  has  disappeared. 

There  exists  an  interesting  relation  between  I"  and  f(w),  which  can  be  found 

n 

as  follows.  In  Eq.  (3.9)  we  substitute  exp(iLgt)Og  for  a,,  and  integrate  the 
result  over  t.  With  Eq.  (A3)  we  then  obtain 


r  iLsT 

rM  -  dt  r(t)  e  a  .  (A5) 

■'o 


as  an  operator  identity.  Then  we  notice  that  Eq.  (3.8)  can  be  inverted  as 


r(0 


1  C  .  -  iwT 

J.i“  ' 


f  (w) 


(A6) 


for  t  >  0.  Substitution  into  Eq.  (AS)  and  performing  the  t  integration  then 
leads  to 


r 


M 


dw  F(w) 


(A7) 


With  the  projectors  of  Section  VII  we  can  write  L^-w  as 


». „  1<ViVa ..  v. 
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(Ls’w)os  •  l  U.b’“)Vs,b  • 


(A8) 


«b 


and  taking  matrix  elements  of  both  sides  gives 


<*,a|((Lg-u»)0g)|b,8>  *  (A-h-w)<a,o|o(;Jb,3> 


ab 


(A9) 


This  shows  that  the  Liouvillian  L^-w  is  diagonal  with  respect  to  the  eigenstates 


of  H-,  and  that  its  matrix  elements  are  A  . -w.  Therefore,  its  inverse  1 / ( L_ ) 

o  MD  b 


has  corresponding  matrix  elements  1/(6^*“)*  which  gives  the  expansion  of 


l/(Lg-uj)  in  projectors  as 


V"  °s 


/  A  -ui  a 

ab 


°SPb  • 


(A10) 


ab 


If  we  insert  this  into  Eq.  ( A 7 )  and  remember  the  general  property 


h  fdw'  ^7  *(u/)  ■  *(u,)  • 

J  -  m 


(  A 1  1  ) 


for  any  Laplace  transform  g(u>),  we  finally  obtain 


fH°S  •  l  f<*.b)(F.«sV 


(A12) 


ab 


Another  way  to  derive  Eq.  ( A 1 2 )  is  by  substituting  the  expansion  (7.6)  for 
exp(iLgt)  into  Eq.  (A5)  and  performing  the  t  integration.  Equation  (A12)  reveals 
that  the  relaxation  operator  in  the  Markov  approximation  effectively  filters  out 
these  w-values  in  f(ui)  which  are  in  exact  resonance  with  the  system  frequencies. 
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